STATE  DOCUMENTS  COLLECTION 

FEB  1  2  2001  ' 

MONTANA  STATE  LIBRARY 

1515  E.  6th  AVE. 
HELENA,  MONTANA  59620 


c&foijtaqa  TisI^, 
J Wildlife  6&<Par*<s 


INVESTIGATION  OF  TRACE  METALS  IN 

AQUATIC  INSECTS  AND  BED  SEDIMENTS  OF  THE 

LOWER  CLARK  FORK  RIVER, 

1996-1997 


by 


Don  Skaar  and  Kurt  Hill 

Montana  Fish,  Wildlife  and  Parks 

1420  East  6th  Avenue 

Helena,  MT  59620 

for  the 

U.S.  Environmental  Protection  Agency 

and 
Montana  Department  of  Environmental  Quality 

March  1999 


MONTANA  STATE  LIBRARY 


..  ...1  tHRI  Inl  llli  IHl — 
3  0864  0015    6686  1 


Investigation  of  trace  metals  in  aquatic  insects  and  bed  sediments  of  the 
lower  Clark  Fork  River,  1996-1997. 

by 

Don  Skaar  and  Kurt  Hill 

Montana  Fish,  Wildlife  and  Parks 

1420  East  6th  Avenue 

Helena,  Montana  §9620 

for  the 

U.S.  Environmental  Protection  Agency 

and 

Montana  Department  of  Environmental  Quality 


March  1999 


Table  of  Contents 

Executive  Summary 3 

Purpose  and  Scope 4 

Methods  and  Materials 4 

Sample  Design  and  Collection 4 

Sample  Analysis  and  Treatment  of  Data 5 

Quality  Assurance  Testing 6 

Results  and  Discussion 6 

Metals  in  fine-grained  sediments 6 

Spatial  Patterns , 9 

Temporal  Patterns 9 

Metals  in  benthic  insects 22 

Spatial  patterns 22 

Temporal  patterns 33 

Community  Composition  of  benthic  invertebrates 46 

Comparison  with  USGS  Data 46 

Acknowledgments 47 

References 49 


Executive  Summary 

Fine-grained  bed  sediments  and  benthic  macroinvertebrates  were  collected  from  three  sites  on 
the  Clark  Fork  River  (East  Missoula,  Missoula,  Alberton)  on  six  occasions  (July  ,  September  and 
November  1996;  March  and  August  1997).  Sediment  and  insect  samples  were  analyzed  for  trace 
elements  believed  to  be  contaminants  at  this  site.  Insect  species  for  which  analyses  were  conducted 
included  Hydropsyche  cockerelli,  H.  occidentalis,  Pteronarcys  californica,  Claassenia  sabulosa,  and 
Arctopsyche  grandis. 

The  mean  concentrations  of  copper  and  zinc  in  fine-grained  sediments  were  highest  at  the 
Missoula  site.  The  concentrations  of  both  metals  at  Missoula  were  not  significantly  different  from 
concentrations  at  East  Missoula,  while  both  sites  were  significantly  different  (and  greater)  than 
concentrations  at  Alberton.  The  East  Missoula  site  showed  no  temporal  trends  for  either  metal, 
while  at  both  the  Missoula  and  Alberton  sites,  concentrations  of  both  copper  and  zinc  were 
significantly  different  (and  lower)  in  July  1996  than  those  measured  in  August  1997. 

The  East  Missoula  site  had  significantly  greater  mean  copper  and  zinc  concentrations  than  the 
Missoula  and  Alberton  sites  for  all  insect  species,  with  the  exception  of  Claassenia,  which  was 
similar  at  all  sites.  In  no  case  was  the  mean  copper  or  zinc  concentration  at  Missoula  significantly 
different  than  Alberton.  Copper  concentrations  in  Hydropsyche  cockerelli  were  higher  in  similarly- 
aged  animals  in  August  1 997  than  in  July  1 996  at  East  Missoula  and  Alberton,  but  slightly  lower  at 
Missoula.  Copper  concentrations  in  Arctopsyche  were  higher  in  similarly-aged  animals  in  August 
1 997  than  in  July  1 996  for  both  East  Missoula  and  Missoula;  this  was  opposite  the  effect  seen  at 
Alberton.  For  Pteronarcys,  copper  concentrations  in  first-  and  second-year  animals  increased  from 
March  to  August  1997,  a  time  period  spanning  the  spring  runoff  flows. 

The  seasonal  pattern  of  benthic  invertebrate  abundance  was  determined  for  the  Missoula  site  and 
results  showed  this  pattern  to  be  similar  to  that  expected  in  temperate-region  streams.  Total  number 
of  organisms  averaged  257.0  per  sample  in  July  1996,  rose  to  423.3  in  September  and  reached  a  peak 
of  809.7  in  November.  Numbers  dropped  substantially  in  January  (176.7)  before  rising  to  399.7  in 
August.  The  precipitous  drop  between  November  and  January  may  have  been  due  to  the  severe  ice 
conditions  and  scour  occurring  at  the  site  during  this  time.  The  seasonal  pattern  of  abundance  of 
individual  insect  taxa  was  also  similar  to  what  would  be  expected. 
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Purpose  and  Scope 

Large  quantities  of  ice  formed  on  the  upper  Clark  Fork  and  Blackfoot  Rivers  in  the  winter  of 
1995-1996.  A  thaw  occurred  in  early  February  1996  and  ice  began  to  break  up  and  move  down  the 
Clark  Fork  River.  Montana  Power  Company  personnel  at  the  Milltown  Dam  quickly  removed  the 
flashboards  on  the  spillway  in  order  to  prevent  damage  to  the  superstructure  of  the  dam.  This 
lowered  the  water  level,  causing  significant  quantities  of  sediment  to  be  washed  from  the  reservoir. 
In  addition,  ice  scour  along  the  bed  and  banks  of  the  Clark  Fork  River  immediately  above  the 
reservoir  contributed  to  the  sediment  load.  Trace  metals  in  the  river  water  below  the  dam  were 
measured  after  this  event  and  found  to  be  quite  high,  particularly  in  the  total  recoverable  fraction. 
On  February  9,  the  U.S.  Geological  Survey  (USGS)  measured  copper  in  the  water  at  concentrations 
of  400  ug/L  (total  recoverable)  and  1 1  ug/L  (dissolved)(USEPA  1998).  Montana  Fish  Wildlife  and 
Parks  (MFWP)  was  concerned  about  the  effects  of  these  metals  on  the  aquatic  ecosystem  of  the  river 
below  Milltown,  and  initiated  a  one  year  study.  The  first  sampling  session  in  July  1 996  was  intended 
to  determine  if  the  metals  mobilized  by  this  ice  event  had  a  lingering  effect  on  the  insects  and 
sediment  of  the  river.  Subsequent  sampling  between  July  1996  and  August  1997  was  intended  to 
establish  reference  levels  of  contamination  from  which  the  effects  of  future  hydrological  events 
could  be  evaluated. 


Methods  and  Materials 

Sample  Design  and  Collection.  The  sampling  design  was  based  upon  the  objective  of  quantifying 
temporal  and  spatial  trends  in  trace  metal  concentrations  in  biota  and  sediments  of  the  Clark  Fork 
River  below  Milltown  Dam.  Three  sites  were  chosen  for  this  purpose:  1)  the  East  Missoula  site  near 
the  Sha-Ron  Fishing  Access  Site  in  East  Missoula  (T13N,  R19W,  secl3).  Exact  location  of  sample 
collections  varied,  but  all  were  within  the  area  from  the  boat  ramp  upstream  about  600  feet  to  the 
pipeline  crossing.  This  site  is  the  same  location  as  the  "Clark  Fork  above  Missoula"  site  of  the 
USGS  (Hornberger  et  al.  1997)  and  "Station  15.5"  of  the  Montana  Department  of  Environmental 
Quality  (McGuire  1 997);  2)  the  Missoula  site  was  near  the  Kelly  Island  Fishing  Access  site  in 
Missoula  (T13N,  R20W,sec20).  Exact  location  of  samples  varied—invertebrates  were  taken 
anywhere  from  the  mouth  of  the  Bitterroot  River  to  about  1000  feet  above  the  access  site.  Sediments 
were  sampled  at  the  access  site  in  July  1 996,  but  thereafter  were  collected  upstream  between  the 
Reserve  Street  Bridge  and  the  Missoula  Wastewater  Treatment  plant;  3)  the  Alberton  site  was 
located  downstream  of  the  town  of  Alberton  (T15N,  R23W,  sec32).  This  site  is  the  same  as  the 
"Clark  Fork  River  below  Missoula"  station  of  the  USGS  (Hornberger  et  al  1997).  The  bulk  of 
invertebrates  were  collected  from  the  same  riffle  area  used  by  the  USGS,  but  occasionally  a  stretch 
of  river  a  mile  above  and  below  this  station  was  used.  Sediments  were  sampled  from  this  station  in 
July  and  September  1996,  but  thereafter  were  sampled  about  two  miles  upstream  at  the  Natural  Pier 
Fishing  Access  Site. 

Bed  sediments  and  benthic  macroinvertebrates  were  collected  for  metals  analysis  using  the 
methods  and  procedures  of  Dodge  et  al.  (1997).  Three  bed  sediment  samples  were  collected  from 


each  site.  Sediments  were  collected  from  the  surface  of  depositional  areas  along  the  river  margin, 
usually  in  the  slow  areas  behind  rocks.  Acid-washed  polypropylene  scoops  were  used  to  collect  the 
sediment.  Sediments  that  appeared  to  be  in  a  reduced  state  (darker  color)  were  not  sampled. 
Sediments  were  passed  through  a  63  urn  filter,  using  river  water  as  a  carrier,  and  stored  in  500  mL 
plastic  bottles. 

Benthic  macroinvertebrates  were  usually  collected  on  the  same  day  as  the  sediments.  Five 
benthic  insect  species  were  targeted  for  metals  analysis:  three  caddisflies  {Hydropsyche  cocker elli, 
Hydropsyche  occidentalis,  Arctopsyche  grandis)  and  two  stoneflies  (Pteronarcys  californica,  and 
Claassenia  sabulosa).  These  taxa  were  selectively  removed  from  kick-net  samples  and  placed  in 
acid-rinsed  plastic  containers  filled  with  river  water.  The  taxa  were  held  together  in  the  containers; 
separate  containers  were  not  provided  for  individual  species.  The  containers  were  placed  on  ice  for 
6-10  hours  to  allow  the  insects  to  depurate  the  contents  of  their  digestive  tracts.  After  this  time,  the 
water  was  drained  and  the  insects  immediately  frozen.  To  determine  community  composition, 
benthic  macroinvertebrates  were  collected  with  a  modified  Hess  sampler  (0. 1  square  meter  diameter, 
1  mm  mesh),  using  methods  described  by  McGuire  (1997).  Three  replicate  samples  were  collected 
from  each  site,  selecting  the  most  heterogenous  habitat  in  riffles.  Samples  were  placed  in  1 -quart 
Mason  jars  and  preserved  with  ethanol.  Samples  were  collected  at  all  sites  on  five  occasions 
between  July  1996  and  August  1997.  The  only  exception  to  this  was  January  1997,  when  ice 
conditions  precluded  sampling  at  all  sites  but  Missoula. 

Sample  Analysis  and  Treatment  of  Data.  For  metals  analysis,  insects  samples  were  thawed  and 
individuals  were  identified  and  sorted  by  species  and  size  groups.  There  were  usually  50-100 
individuals  per  size  group  for  the  caddisfly  species,  while  1 0-20  individuals  was  more  typical  for  the 
stoneflies.  For  the  largest  size  groups  of  stoneflies,  there  were  occasions  where  as  few  as  five 
individuals  were  used.  For  the  largest  Claassenia  (greater  than  25  mm  in  length)  there  was  probably 
no  uncertainty  with  regard  to  correct  identification,  but  with  smaller  animals  there  was  probably 
occasional  misidentification  of  a  few  individuals  which  may  have  been  any  one  of  the  following 
genera:  Isogenoides,  Skwala,  Calineuria.  Ten  individuals  from  each  size  group  were  measured  for 
total  length.  Each  insect  sample  was  rinsed  with  distilled  water,  transferred  to  acid-washed  plastic 
centrifuge  vials  and  placed  in  an  oven  at  60°C  until  dry.  These  samples  were  then  ground  with  a 
mortar  and  pestle  and  prepared  for  metals  analysis  by  microwave-digestion  for  approximately  6 
minutes  in  a  solution  of  4  mL  nitric  acid  and  1  mL  hydrochloric  acid.  The  digestate  was  then 
brought  to  a  50  mL  volume  with  distilled  water.  Samples  were  then  analyzed  by  graphite  furnace 
AAS  (Atomic  Absorption  Spectroscopy)  or  ICP  (Inductively  Coupled  Plasma)  analysis.  The 
procedure  for  metals  analysis  of  sediments  was  similar  to  insects,  except  that  sediments  were  dried 
in  the  same  plastic  vials  in  which  they  had  been  collected. 

Statistical  calculations  and  tests  were  performed  by  using  SPSS  software  (SPSS,  1993). 
Determinations  of  statistical  significance  were  based  on  PO.05  levels. 


Quality  Assurance  Testing 

Procedural  blanks  for  bed  sediment  and  insect  samples  contained  the  same  reagents  (nitric  and 
hydrochloric  acid,  distilled  water)  that  were  used  for  sample  preparation.  Zinc  was  the  only  element 
for  which  contamination  was  detected  (Table  1).  One  of  these  contaminated  samples  had  a 
concentration  of  0.006  ug/mL,  which  is  very  near  to  the  detection  limit  (0.005  ug/mL),  while  the 
other  sample  was  much  higher  (0.032  ug/mL).  This  highly  contaminated  blank  actually  contained 
0.32  ug  zinc  (0.032  ug/mL  X  10  mL),  which  is  about  4%  of  the  amount  of  zinc  in  the  insect  samples 
in  this  study  (80  ug/g  zinc  X  0.10  g  dry  weight  =  8  ug  zinc)  or  1%  of  the  amount  of  zinc  in  the 
sediment  samples  (120  ug/g  zinc  X  0.25  g  dry  weight  =  30  ug  zinc).  These  errors  are  considered 
minor  and  no  adjustment  was  made  to  any  of  the  insect  or  bed-sediment  values. 

Standard  reference  materials  used  in  this  study  were  bovine  liver  (as  a  surrogate  for  insect  tissue) 
and  sediment.  Mean  recovery  efficiency  for  bovine  liver  was  within  90%  of  the  certified 
concentration  for  copper  and  cadmium,  but  was  low  for  zinc  (68%)  and  high  for  lead  (198%)  and 
arsenic  (206%)(Table  2).  Recovery  for  sediments  was  much  better;  mean  recovery  efficiency  ranged 
from  97%  for  arsenic  to  105%  for  zinc.  The  95%  confidence  intervals  for  recovery  efficiency  were 
always  greater  than  the  error  considered  acceptable.  The  certified  acceptance  limits  were  about  90 
and  110%  of  the  mean  value  for  arsenic,  cadmium,  lead  and  zinc,  and  60  and  140%  for  copper.  In 
comparison,  the  confidence  levels  that  were  achieved  were  about  80  and  120%  of  the  mean  value 
for  arsenic,  cadmium,  lead  and  zinc,  and  55  and  145%  for  copper. 


Results  and  Discussion 

Metals  in  fine-grained  sediments.  An  understanding  of  the  trends  in  sediment  metals  must  include 
an  examination  of  the  hydrologic  events  that  transported  sediments  to  the  sites,  and  the  extent  to 
which  the  sampled  sediments  were  reflective  of  those  events.  The  discharges  occurring  at  our  three 
sites  are  fairly  well  portrayed  by  measurements  taken  at  the  USGS  gage  stations  above  Missoula 
(near  East  Missoula)  and  below  Missoula  (below  the  confluence  of  the  Bitterroot  River).  The 
discharges  measured  near  East  Missoula  should  be  very  similar  to  those  occurring  at  our  East 
Missoula  and  Missoula  sites,  while  the  discharges  from  the  station  below  Missoula  are  somewhat 
less  than  those  occurring  at  our  Alberton  site.  At  both  stations,  the  peak  flows  from  the  February 
1996  ice-out  occurred  on  the  10th,  when  flows  averaged  12,900  cfs  above  Missoula  and  20,300  cfs 
below  Missoula.  Over  3 1 ,000  tons  of  sediment  were  discharged  from  Milltown  Dam  during  the  four 
peak  days  of  the  February  event  (USEPA  1998).  These  sediments  were  probably  transported  farther 
downstream  as  a  result  of  the  higher  flows  during  spring  runoff  in  1996.  These  flows  peaked  at 
18,300  cfs  on  June  10  above  Missoula  and  at  37,400  cfs  on  May  19  below  Missoula.  After  the  spring 
runoff  in  1 996,  the  hydrograph  followed  a  typical  pattern  of  flows  dropping  through  the  fall  and 
winter.  High  flows  were  not  seen  on  the  river  again  until  spring  1997.  The  1997  flows  were  much 
higher  than  in  1996  and  peaked  on  May  18  at  both  stations  when  discharge  got  to  26,400  cfs  above 
Missoula  and  54,100  cfs  below  Missoula. 


Table  1.  Metal  analyses  of  procedural  blanks  for  insect  and  sediment  samples.  Changes  in  detection 
limits  for  a  specific  element  reflect  a  change  in  instrument  used  (ICP  vs.  AAS). 

Sample  Period 

Dilution 
Ratio 

Metal  concentration  (ug/mL) 

Arsenic 

Cadmium 

Copper 

Lead 

Zinc 

July  1996 

1:9 

<0.001 

O.001 

<0.005 

O.001 

<0.005 

July  1996 

1:9 

O.001 

<0.001 

<0.005 

<0.05 

O.005 

July  1996 

1:9 

<0.001 

<0.001 

<0.005 

<0.001 

<0.005 

September  1996 

1:9 

<0.001 

<0.001 

<0.005 

<0.001 

<0.005 

September  1996 

1:9 

O.001 

O.001 

<0.005 

<0.05 

<0.005 

September  1996 

1:9 

O.001 

O.001 

<0.005 

<0.05 

<0.005 

November  1996 

1:9 

O.001 

<0.001 

<0.005 

<0.05 

O.005 

November  1996 

1:9 

O.001 

<0.001 

O.005 

<0.05 

0.032 

March  1997 

1:9 

<0.001 

<0.002 

<0.005 

<0.001 

0.006 

March  1997 

1:9 

<0.001 

<0.001 

<0.005 

<0.05 

<0.005 

August  1997 

1:9 

O.001 

O.001 

<0.005 

<0.05 

<0.005 

August  1997 

1:9 

<0.001 

<0.001 

<0.005 

<0.05 

<0.005 

Table  2.  Recovery  efficiency  for  metal  analysis  of  standard  reference  materials  (SRM)  for 
insects  and  sediment. 

Metal 

Number  of 
measurements 

Certified 

concentration 

(ug/g) 

SRM  recovery  (ug/g) 
Mean        95%  CI 

Acceptance 
limits  for  SRM 
(ug/g) 

Insect  (NBS  #1577,  Bovine  Liver) 

Arsenic 

5 

0.050* 

0.103 

0.025-0.181 

..** 

Cadmium 

5 

0.50 

0.47 

0.43-0.51 

- 

Copper 

5 

160 

144 

132.8-149.6 

~ 

Lead 

4 

0.13 

0.257 

0.217-0.297 

-- 

Zinc 

4 

127 

86.7 

74.7-98.7 

- 

Sediment  (USGS  #5) 

Arsenic 

6 

216 

210 

166-256 

196-236 

Cadmium 

6 

158 

162 

128-196 

142-174 

Copper 

4 

58.9 

66.6 

34.1-97.1 

31.3-82.5 

Lead 

7 

309 

324 

258-390 

291-327 

Zinc 

7 

598 

588 

500-676 

551-645 

""Uncertified  concentration 
**Ranges  not  provided  with  material 


All  sediments  for  this  study  were  sampled  at  discharges  between  1,500  and  7,500  cfs,  which 
correspond  to  water  surface  elevations  5-10  feet  lower  than  those  during  spring  runoff  and  3-8  feet 
lower  than  during  the  1996  ice-out  (Table  3).  Therefore,  all  of  the  sediment  samples  were  clearly 
collected  from  within  the  portion  of  the  channel  that  had  been  influenced  by  all  of  these  events.  The 
fact  that  spring  runoff  flows  in  1 996  and  1 997  were  much  higher  than  during  the  ice-out  event 
strongly  suggests  that  the  ice-out  sediments  were  being  transported  through  and  downstream  from 
our  sites  during  the  course  of  our  study.  However,  we  have  no  way  to  know  when  or  where  we  may 
have  sampled  the  ice-out  sediments. 

Spatial  Patterns.  The  mean  concentrations  of  most  metals  in  the  sediments  (arsenic,  copper,  iron, 
manganese,  nickel,  lead  and  zinc)  were  highest  at  the  Missoula  site.  Mean  concentrations  of  silver 
and  cadmium  were  highest  at  East  Missoula,  while  concentrations  of  chromium  were  highest  at 
Alberton  (Tables  4  and  5).  Concentrations  of  copper  and  zinc  at  Missoula  averaged  270.5  and  556.4 
ug/g,  respectively,  only  slightly  higher  than  the  corresponding  concentrations  at  East  Missoula  (247.8 
and  531.9  ug/g),  but  much  higher  than  at  Alberton  (129.7  and  274.8  ug/g).  The  data  for  these  two 
metals  were  subjected  to  statistical  analysis,  and  results  showed  that  East  Missoula  and  Missoula 
were  not  significantly  different  in  terms  of  mean  copper  or  zinc  concentrations  (Table  6),  while  both 
of  these  sites  were  significantly  different  (and  greater)  than  concentrations  at  Alberton. 

Temporal  Patterns.  There  were  distinct  differences  between  the  sites  in  terms  of  temporal 
changes  in  metals  in  the  sediments.  One-way  ANOVA  was  used  to  test  the  null  hypothesis  that  the 
metal  concentrations  on  all  sampling  dates  had  equal  means.  At  East  Missoula,  no  significant 
differences  were  seen  for  any  of  the  metals  tested  (Figures  1-4,  Table  7).  In  contrast,  the  null 
hypothesis  was  rejected  for  arsenic,  cadmium,  copper,  lead,  and  zinc  at  both  the  Missoula  and 
Alberton  sites.  Measurements  of  copper  and  zinc  at  these  two  sites  were  subjected  to  further 
statistical  analysis.  Bonferroni  t-test  comparisons  were  done  for  the  mean  concentration  of  copper 
and  zinc  on  each  sampling  date.  Where  significant  differences  were  found,  almost  all  showed  an 
increase  in  concentration  over  time.  At  Missoula,  all  of  the  differences  (10  of  10  tests)  showed  an 
increase  over  time,  while  at  Alberton,  1 1  of  14  significant  tests  showed  an  increase  over  time  (Tables 
8  and  9).  At  both  sites,  copper  and  zinc  levels  measured  in  July  and  September  1996  were 
significantly  different  (and  lower)  than  those  measured  in  August  1 997.  These  findings  are  qualified 
by  the  fact  that  the  sampling  location  for  sediments  was  changed  for  both  sites  during  the  study  (after 
July  1 996  at  Missoula  and  after  September  1 996  at  Alberton),  and  it  is  possible  that  there  were 
differences  in  metal  concentrations  between  the  two  locations.  Nonetheless,  the  September  1 996 
levels  of  copper  and  zinc  were  still  significantly  lower  at  Missoula  than  the  August  1997  samples; 
similarly  the  November  1996  copper  and  zinc  levels  at  Alberton  were  significantly  lower  than  in 
August  1997. 

One  possible  explanation  for  the  increase  in  copper  and  zinc  levels  over  time  at  Missoula  and 
Alberton,  but  not  East  Missoula,  is  that  it  was  a  reflection  of  the  movement  of  the  sediments  that  had 
washed  out  of  Milltown  Dam  in  February  1996.  This  explanation  assumes  the  sediment  from  the 
reservoir  had  moved  past  the  East  Missoula  site  by  the  time  this  study  began  in  July  1 996  but  was 


Table  3.  Water  surface  elevations  (ft)  at  two  Clark  Fork 

River  Stations  during  1996  and  1997.  Elevations  are 
relative  to  datum  at  each  site,  and  represent  mean 

daily  elevations  unless  noted  otherwise.  (USGS 

data). 

Clark  Fork                       Clark  Fork 
Date               above  Missoula             below  Missoula 

2-10-96 

9.04  (max.  elev.) 

7.54  (max.  elev.) 

5-19-96 

- 

9.20  (max.  elev.) 

6-10-96 

9.94  (max.  elev.) 

- 

7-29-96 

3.40 

~ 

7-31-96 

~ 

2.48 

9-3-96 

2.80 

~ 

9-4-96 

- 

1.70 

11-5-96 

2.92 

- 

11-6-96 

- 

1.78 

3-25-96 

4.60 

- 

3-26-96 

~ 

4.00 

5-18-97 

12.59  (max.  elev.) 

12.18  (max.  elev.) 

8-12-97 

3.45 

2.63 
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Table  6.    Comparisons  of  mean  copper  and  zinc  concentrations  in  fine-grained  sediments  at 

the  three  sites.  Levels  of  significance  (P)  are  based  on  Bonferroni  t-test  comparisons. 

Dependent 

variable                      Site  A                       Site  B 

Mean  Difference 
(A-B) 

Sig. 

Copper 

East 
Missoula 

Missoula 

-22.7 

.856 

East 
Missoula 

Alberton 

118.1* 

.000 

Missoula 

Alberton 

140.8* 

.000 

Zinc 

East 
Missoula 

Missoula 

-24.5 

1.000 

East 
Missoula 

Alberton 

257.1* 

.000 

Missoula 

Alberton 

281.6* 

.000 

*The  mean  difference  is  significant  at  the  P=0.05  level. 


14 


P>     400 


c 
g 

15 


c 

0) 

u 

c 
o 
u 

d) 

Q. 
Q. 
O 

o 

E 

=3 

E 


E 

3 

E 
x 

(0 

E 


(0 


300^ 


2004 


100- 


SITE 


~~  East  Missoula 


Missoula 


Alberton 


7/96        9/96       11/96       1/97        3/97        8/97 


Figure  1 .  Copper  concentrations  (ug/g  dry  weight)  in  fine-grained  bed 
sediments  of  the  Clark  Fork  River  at  three  sites. 
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Figure  2.  Zinc  concentrations  (ug/g  dry  weight)  in  fine-grained  bed 
sediments  of  the  Clark  Fork  River  at  three  sites. 
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Figure  3.  Cadmium  (ug/g  dry  wt)  in  fine-grained  bed  sediments.  Negative 
value  indicates  that  concentration  was  below  detection  limit. 
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Figure  4.  Lead  concentrations  (ug/g  dry  wt)  in  fine-grained  bed 
sediments  of  the  Clark  Fork  River  at  three  sites. 


Table  7.  Results  of  one-way  ANOVA  tests  of  the  differences  between  mean  metal 
concentrations  in  fine-grained  sediments  at  each  site  over  time.  In  the 
case  where  homogeneity  of  variance  test  failed,  significance  level  is  based 
on  the  Kruskal  Wallis  (K-W)  nonparametric  test. 

East  Missoula 

Missoula 

Alberton 

F 

Sig. 

F 

Sig. 

F 

Sig. 

Arsenic 

0.18          0.942 

4.46             0.018* 

80.3             0.000* 

Cadmium 

0.30          0.870 

59.1             0.000* 

K-W           0.014* 

Copper 

0.76          0.574 

10.1             0.001* 

56.4            0.000* 

Manganese 

0.59          0.676 

K-W           0.834 

0.67            0.631 

Lead 

1.25           0.351 

5.32            0.010* 

3.99            0.035* 

Zinc 

K-W         0.271 

15.6            0.000* 

42.2            0.000* 

♦Significant  at  P<0.05  level 
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still  moving  into  and  through  the  other  sites.  However,  it  is  also  possible  that  these  trends  were  due 
to  or  influenced  by  unique  or  different  conditions  at  each  site  which  might  affect  the  metal 
concentrations.  Boggs  (1994)  suggested  that  such  conditions  might  include  differences  in 
streamflow,  sediment  deposition,  sediment  size  fractions,  organic  matter  content,  metal-oxide 
content,  and  proportions  of  oxidized  to  reduced  sediments.  A  further  concern  is  that  we  may  have 
committed  Type  I  or  Type  II  errors  in  our  statistical  analysis;  the  error  being  due  to  the  inherent 
variability  of  sediment  metal  concentrations  or  perhaps  an  inappropriate  temporal  scale  within  which 
we  sampled  (2-4  month  intervals). 

The  high  flows  during  spring  runoff  in  1 997  apparently  resulted  in  a  much  larger  loading  of 
metals  to  the  Clark  Fork  River  than  had  been  recorded  during  the  previous  12  years  (Lambing  1998). 
However,  based  on  samples  taken  for  this  study,  there  was  no  clear  effect  on  concentrations  of 
copper  and  zinc  in  sediments.  As  discussed  above,  there  were  no  significant  changes  over  time  at 
East  Missoula.  Even  at  Missoula  and  Alberton,  where  the  ANOVA  test  revealed  significant  changes, 
there  was  no  difference  in  mean  zinc  concentrations  at  either  site  when  comparing  pre-runoff 
conditions  (March  1997  samples)  to  post-runoff  conditions  (August  1997).  For  copper,  no  difference 
was  seen  at  the  Missoula  site  for  the  two  time  periods,  while  at  Alberton,  copper  actually  decreased 
from  March  to  August  1997. 

Metals  in  benthic  insects.  Although  spatial  and  temporal  trends  were  seen  with  all  metals,  the  only 
metals  subjected  to  statistical  analysis  were  copper  and  zinc. 

Spatial  Patterns.  The  Hydropsyche  species  (H.  cockerelli  and  H.  occidentalis)  displayed  the 
most  consistent  drop  in  copper  concentrations  from  East  Missoula  to  Missoula  and  then  to  Alberton. 
For  all  dates  combined,  the  mean  copper  concentration  was  55.0  and  52.9  ug/g  at  East  Missoula, 
decreasing  to  40.5  and  35.6  ug/g  at  Missoula,  and  29.7  and  28.3  ug/g  at  Alberton  for  H.  cockerelli 
and  H.  occidentalis,  respectively  (Tables  10,  11  and  12).  For  Pteronarcys,  mean  copper 
concentrations  were  highest  at  East  Missoula  (71.0  ug/g)  but  similar  at  Missoula  (52.4  ug/g)  and 
Alberton  (53.0  ug/g)(Table  13).  Claassenia  showed  the  least  variation  in  copper  between  sites, 
ranging  from  a  high  of  48.4  ug/g  at  East  Missoula  to  a  low  of  44.3  ug/g  at  Alberton  (Table  14). 
Arctopsyche  had  the  lowest  copper  concentrations  of  all  species,  decreasing  from  a  high  of  40.0  ug/g 
at  East  Missoula  to  a  low  of  1 8.0  ug/g  at  Alberton  (Table  1 5). 

The  spatial  patterns  of  zinc  were  generally  similar  to  those  of  copper.  Claassenia  had  the  highest 
concentrations  of  all  species,  which  were  also  low  in  variability  between  sites  (ranging  from  a  high 
of  179.3  ug/g  at  East  Missoula  to  a  low  of  171 .7  ug/g  at  Missoula).  Pteronarcys  also  had  high  zinc 
levels  at  East  Missoula  (178.1  ug/g),  but  concentrations  decreased  to  their  lowest  at  Missoula  (153.9 
ug/g).  The  Hydropsyche  species  were  quite  similar:  both  had  highest  zinc  levels  at  East  Missoula 
(172.8  and  173.3  ug/g  for  H.  occidentalis  and  H.  cockerelli,  respectively)  and  lowest  levels  at 
Alberton  (136.0  and  121.2  ug/g).  Arctopsyche  had  the  lowest  zinc  concentrations,  but  still  showed 
a  decrease  in  a  downstream  direction  (148.0  ug/g  at  East  Missoula,  121.3  ug/g  at  Missoula  and  119.9 
ug/g  at  Alberton). 
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Statistical  analysis  of  these  data  with  one-way  ANOVA  showed  that  for  all  species  but 
Claassenia  we  can  reject  the  null  hypothesis  (at  the  P=0.05  level)  that  the  mean  copper 
concentrations  are  the  same  at  the  three  sites  (Table  16).  For  zinc,  the  null  hypothesis  was  rejected 
in  Pteronarcys,  H.  Cockerelli,  and  Arctopsyche,  but  not  H.  occidentalis  and  Claassenia.  In  all  cases 
where  the  null  hypothesis  was  rejected,  it  was  because  of  differences  between  the  high  mean  values 
at  the  East  Missoula  site  and  lower  values  at  one  of  the  other  sites.  Pairwise  comparisons  showed 
that  the  East  Missoula  site  had  a  significantly  greater  mean  copper  value  than  Missoula  for 
Pteronarcys,  both  Hydropsy 'c he  species,  and  Arctopsyche  (Table  17).  Pteronarcys,  H.  cockerelli  and 
Arctopsyche  also  had  significantly  higher  mean  concentrations  of  zinc  at  East  Missoula  than  at 
Missoula.  Copper  concentrations  were  also  significantly  higher  at  East  Missoula  than  Alberton  for 
both  Hydropsyche  species  and  Arctopsyche;  for  zinc,  only  H.  cockerelli  was  significantly  higher  at 
East  Missoula  than  Alberton.  In  no  case  was  the  mean  copper  or  zinc  concentrations  at  Missoula 
significantly  different  than  Alberton. 

Temporal  Patterns.  Insects  were  sampled  five  times  at  East  Missoula  and  Alberton,  and  six 
times  at  Missoula  during  the  thirteen  months  of  the  study.  This  was  done  in  order  to  identify 
changes  in  metal  concentrations  in  insects  that  could  result  from  significant  hydrologic  events  in  the 
drainage  or  changes  in  Milltown  Dam  operation.  This  type  of  evaluation  must  consider  the  periods 
of  emergence  of  the  insects  being  sampled  and  account  for  the  different  cohorts  for  species  living 
longer  than  one  year  (e.g.  three  years  for  Pteronarcys).  Consequently,  we  have  restricted  our 
analysis  of  temporal  trends  to  those  situations  where  the  age  of  the  insects  is  clearly  understood.  We 
have  used  length-frequency  analysis  as  the  technique  for  estimating  the  age  of  animals.  For 
example,  H.  cockerelli  larva  at  the  Missoula  site  were  close  to  emerging  on  March  27,  1997  when 
lengths  ranged  from  13.0-19.5  mm.  When  this  species  was  sampled  five  months  later  on  August  12, 
1997,  the  lengths  ranged  from  10.0-13.0  mm,  which  we  assumed  to  be  from  a  new  generation  of 
animals.  Also,  we  have  not  included  Claassenia  in  this  analysis  because  it  may  have  some  ability 
to  regulate  copper  (Cain  et  al.  1992)  and  therefore  seems  to  be  a  poor  barometer  of  changes  in  the 
copper  levels  in  the  physical  media  (water,  sediment). 

Pteronarcys  was  the  only  species  for  which  we  were  able  to  sample  individuals  that  had  lived 
through  the  February  1 996  ice-out  event.  These  individuals  were  first  sampled  in  July  1 996  as  2nd- 
and  3rd-year  animals.  A  comparison  of  the  metal  concentrations  in  these  animals  with  younger 
animals  born  during  or  after  high  water  in  1 996  can  serve  to  isolate  some  of  the  effects  of  the 
February  event.  Results  show  that  the  lst-year  animals,  sampled  in  September  or  November  1996, 
had  similar  or  higher  levels  of  both  copper  or  zinc  (Table  18).  This  suggests  that  the  older 
Pteronarcys  animals  did  not  have  unexpectedly  high  levels  of  copper  and  zinc,  and  if  there  were  any 
lingering  effects  from  the  February  event,  they  were  not  manifested  as  high  body  burdens  of  metals. 

Pteronarcys  was  also  the  only  species  we  sampled  that  could  be  used  to  evaluate  whether  there 
were  changes  in  copper  and  zinc  concentrations  in  animals  that  lived  through  the  1 997  runoff.  This 
comparison  was  possible  for  a  combination  of  first-  and  second-year  animals  at  Missoula  and  East 
Missoula  (Table  19).  At  East  Missoula,  the  mean  concentrations  of  the  metals  changed  little  from 
March  to  August:  copper  increased  slightly  from  68.8  ug/g  to  80.9  ug/g,  while  zinc  decreased  from 
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Table  16.  Results  of  one-way  ANOVA  tests  of  the  differences  between  mean  copper  and 
zinc  concentrations  at  the  three  sites. 

Pteronarcys 
californica 

Hydropsyche 
cockerelli 

Hydropsyche 
occidentalis 

Arctopsyche 
grandis 

Claassenia 
sabulosa 

F 

Sig. 

F 

Sig. 

F 

Sig. 

F 

Sig. 

F 

Sig. 

Copper 

3.7       0.035 

13.5       0.000 

7.6        0.006 

12.3        0.000 

0.8        0.473 

Zinc 

6.4       0.005 

9.0       0.002 

1.6        0.239 

5.1        0.013 

0.9        0.404 
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Table  18.  Concentration  of  copper  and  zinc  (ug/g  dry  weight)  in 
groups  in  the  Clark  Fork  River,  1996. 

July  1996                                September  1996 
(2nd  yr  or  older  animals)                      ( 1  st  yr  animals) 

different  Pteronarcys  size 

November  1996 
( 1  st  yr  animals) 

Size  range 
(mm) 

Cu            Zn 

Size  range 
(mm) 

Cu        Zn 

Size  range 
(mm) 

Cu        Zn 

East  Missoula 

13.5-23.5 

59.5 

170 

9.5-21.0 

58.1 

184 

7.5-20.5 

81.1 

205 

29.0-37.5 

64.4 

156 

Mean=  62.0             163 

Missoula 

16.5-26.0 

35.1 

146 

10.0-12.5 

48.6 

206 

11.5-17.0 

49.1 

168 

26.5-33.0 

41.1 

138 

34.5-43.0 

51.4 

144 

Mean  =     42.5            143 
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Table  19.  Comparison  of  copper  and  zinc  concentrations  (ug/g  dry  wt)  in  different  size 
groups  of  Pteronarcys  before  and  after  exposure  to  runoff  in  spring  1997. 

March  1997 

August  1997 

Site 

Size  range 
(mm) 

Cu 

Zn 

Size  range 
(mm) 

Cu 

Zn 

East  Missoula 

9.0-13.0 

65 

198 

12.0-17.0 

85.8 

186 

East  Missoula 

21.0-29.0 

72.7 

163 

21.0-26.0 

66.7 

164 

East  Missoula 

32.5-37.0 

90.3 

166 

Mean  = 

68.8 

180.5 

Mean  = 

80.9 

172 

Missoula 

11.0-18.0 

37 

134 

22.0-26.0 

128 

185 

Missoula 

19.0-26.0 

41.1 

140 

30.0-36.0 

102 

193 

Mean  = 

39.0 

137 

Mean  = 

115 

188 
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1 80.5  ug/g  to  1 72  ug/g.  This  was  in  contrast  to  the  trend  at  Missoula  where  copper  concentrations 
increased  markedly  from  39  to  115  ug/g  and  the  zinc  increased  from  137  to  188  ug/g. 

Another  useful  comparison  that  can  be  made  is  of  the  metal  concentrations  in  animals  born 
slightly  before  or  during  the  runoff  periods  of  1996  and  1997.  Any  differences  in  metal 
concentrations  can  probably  be  attributed,  at  least  in  part,  to  differing  amounts  of  exposure  to  metals 
in  water  and  sediment  during  the  runoff  itself.  This  comparison  was  possible  for  Hydropsyche 
cockerelli  and  Arctopsyche  at  all  sites.  Four  of  the  six  comparisons  for  copper  showed  an  increase 
from  1996,  with  an  average  increase  of  32%  for  H.  cockerelli  and  99%  for  Arctopsyche  for  all  three 
sites  combined  (Table  20).  For  zinc  there  was  no  obvious  trend,  as  only  three  of  six  comparisons 
showed  an  increase  over  time.  Furthermore,  the  average  change  in  zinc  concentration  over  time  was 
a  -  17%  for  H.  cockerelli  and  +12%  for  Arctopsyche.  These  results  suggest  that  conditions  were 
more  conducive  for  the  accumulation  of  copper,  but  not  zinc,  during  the  spring  of  1 997  than  the 
spring  of  1996. 

Community  composition  of  benthic  invertebrates.  Three  Hess  samples  were  collected  from  each 
site  in  July,  September  and  November  of  1996  and  again  in  August  of  1997.  Ice  conditions  in 
January  1 997  prevented  us  from  sampling  all  but  the  Missoula  site,  while  rising  water  in  March  1 997 
precluded  sampling  at  all  sites. 

The  Hess  samples  from  the  Missoula  site  were  analyzed  and  summarized  for  this  report.  The 
total  number  of  organisms  averaged  257.0  per  sample  in  July  1996,  rose  to  423.3  in  September  and 
reached  a  peak  of  809.7  in  November.  Numbers  dropped  substantially  in  January  (176.7)  before 
rising  to  399.7  in  August  1997  (Figure  5,  Tables  21-25).  This  pattern  is  generally  consistent  with 
that  described  by  Hynes  (1970)  for  stream-dwelling  insects  in  temperate  regions.  This  pattern  is  one 
in  which  most  insects  reach  their  peak  numbers  in  late  fall-early  winter  as  a  result  of  the  growth  and 
recruitment  of  young  from  reproduction  the  previous  summer.  Throughout  the  winter,  numbers 
decline  as  a  result  of  death,  and  by  spring  additional  reductions  result  from  insect  emergence. 
Minimum  numbers  are  then  typically  reached  in  late  spring.  On  the  Clark  Fork  there  was  a 
precipitous  drop  in  just  two  months  from  November  to  January.  Whether  or  not  this  is  typical  is 
unknown,  but  the  severe  ice  conditions  and  scour  on  the  river  at  the  time  may  have  had  an  effect  on 
the  densities  of  near-shore  dwellers.  Regardless  of  the  cause  for  the  low  numbers  in  January,  the 
populations  subsequently  increased  to  the  point  that  the  August  1997  numbers  (399.7)  were  actually 
somewhat  higher  than  the  July  1996  numbers  (257.0). 

This  seasonal  pattern  in  abundance  described  by  Hynes  (1970)  was  followed  at  the  Missoula  site 
by  Dipterans,  caddisflies  and  stoneflies,  but  not  by  mayflies.  Mayflies  had  a  density  of  174.3  per 
Hess  sample  in  July  1996,  but  then  dropped  to  69.7  in  September  and  back  up  again  to  174.7  in 
November.  No  obvious  explanation  can  be  provided  for  this  pattern.  Caddisflies  were  unusual  in 
that  their  highest  densities  were  reached  in  August  1997,  at  which  time  they  were  more  than  nine 
times  more  plentiful  (259)  than  they  were  in  July  1996  (27.7).  Most  of  the  difference  between  these 
two  samples  was  due  to  the  increase  in  numbers  of  one  species  (H.  cockerelli)  in  August  1997.  On 
that  date,  H.  cockerelli  made  up  67.2%  of  the  total  caddisflies  at  the  site,  whereas  in  July  1996  the 
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Table  20.  Comparison  of  copper  and  zinc  concentrations  (ug/g  dry  wt)  in  insects  born  before 
or  during  spring  runoff. 

July  1996 

August  1997 

Site 

Species 

Size  range 
(mm) 

Cu 

Zn 

Size  range 
(mm) 

Cu 

Zn 

E.  Missoula 

H.  cockerelli 

7.0-11.5 

39.7 

210 

9.0-14.0 

67.8 

136 

Missoula 

H.  cockerelli 

7.0-15.0 

48.6 

180 

10.0-13.0 

46.3 

136 

Alberton 

H.  cockerelli 

6.5-14.0 

28.5 

119 

11.0-14.0 

37.4 

131 

Average  change  1996  to  1997  = 

+32% 

-17% 

E.  Missoula 

Arctopsyche 

10.0-16.0 

19.8 

118 

8.0-11.0 

64.5 

171 

Missoula 

Arctopsyche 

9.0-14.0 

12.6 

110 

12.0-14.5 

24 

123 

Alberton 

Arctopsyche 

6.0-14.0 

20.1 

149 

12.0-17.0 

16.4 

118 

Average  change  1996  to  1997  = 

+99% 

+12% 
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Figure  5.  Number  of  benthic  organisms  per  Hess  sample  from  the 
Missoula  site.  N=4  for  7  &  9/96;  N=3  for  other  dates. 


Table  2 1 .    Enumeration  of  macroinvertebrate  taxa  from  Hess  samples  at  the  Missoula  site,  July  24, 
1996. 

Taxon 

Sample 
1 

Sample 

2 

Sample 

3 

Sample 

4 

SUM 

Mean 

SD 

COLEOPTERA 

0 

4 

5 

0 

9 

2.2 

2.6 

DIPTERA 

11 

26 

150 

0 

187 

46.7 

69.7 

EPHEMEROPTERA 

95 

109 

301 

192 

697 

174.3 

94.7 

PLECOPTERA 

2 

1 

7 

11 

21 

5.2 

4.6 

Claassenia  sabulosa 

1 

0 

5 

3 

9 

2.3 

2.2 

Pteronarcys 
californica 

1 

0 

2 

8 

11 

2.7 

3.6 

TRICHOPTERA 

24 

12 

40 

35 

111 

27.7 

12.4 

Hydropsyche  cockerelli 

2 

1 

6 

2 

11 

2.7 

2.2 

Hydropsyche 
occidentalis 

4 

2 

10 

12 

28 

7.0 

4.8 

Arctopsyche  grandis 

3 

1 

1 

4 

9 

2.2 

1.5 

AMPHIPODA 

0 

0 

1 

0 

1 

0.3 

0.5 

GASTROPODA 

0 

0 

0 

1 

1 

0.3 

0.5 

ANNELIDA 

0 

0 

0 

1 

1 

0.3 

0.5 

TOTAL  ORGANISMS 

132 

152 

504 

240 

1028 

257.0 

171.2 
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Table  22.  Enumeration  of  macroinvertebrate  taxa  from  Hess  samples  at  the  Missoula  site,  September  3,  1996. 

Taxon 

Sample  1 

Sample  2 

Sample  3 

Sample  4 

SUM 

Mean 

SD 

COLEOPTERA 

5 

7 

10 

11 

33 

8.3 

2.8 

DIPTERA 

95 

158 

251 

119 

623 

155.7 

68.6 

EPHEMEROPTERA 

79 

66 

89 

45 

279 

69.7 

19.0 

PLECOPTERA 

17 

17 

15 

24 

73 

18.3 

3.9 

Claassenia  sabulosa 

11 

10 

10 

11 

42 

10.5 

0.6 

Pteronarcys  californica 

6 

7 

5 

10 

28 

7.0 

2.2 

TRICHOPTERA 

90 

173 

268 

152 

683 

170.7 

73.8 

Hydropsyche  cockerelli 

34 

62 

44 

47 

187 

46.7 

11.6 

Hydropsyche  occidentalis 

34 

64 

178 

79 

355 

88.7 

62.4 

Arctopsyche  grandis 

4 

4 

13 

9 

30 

7.5 

4.4 

ANNELIDA 

2 

0 

0 

0 

2 

0.5 

1.0 

TOTAL  ORGANISMS 

288 

421 

633 

351 

1693 

423.3 

150.0 

42 


Table  23.  Enumeration  of  macroinvertebrate  taxa  from  Hess  samples  at  the  Missoula  site,  November  5, 
1996 

Taxon 

Sample  1 

Sample  2 

Sample  3 

SUM 

Mean 

SD 

COLEOPTERA 

6 

3 

20 

29 

9.7 

9.1 

DIPTERA 

232 

415 

317 

964 

321.3 

91.6 

EPHEMEROPTERA 

128 

210 

186 

524 

174.7 

42.1 

PLECOPTERA 

75 

70 

116 

261 

87.0 

25.2 

Claassenia  sabulosa 

22 

12 

26 

60 

20.0 

7.2 

Pteronarcys  californica 

8 

1 

8 

17 

5.7 

4.0 

TRICHOPTERA 

291 

120 

239 

650 

216.7 

87.7 

Hydropsyche  cockerelli 

70 

22 

73 

165 

55.0 

28.6 

Hydropsyche  occidentalis 

171 

64 

111 

349 

115.3 

53.6 

Arctopsyche  grandis 

7 

4 

15 

26 

8.7 

5.7 

GASTROPODA 

0 

0 

1 

1 

0.3 

0.6 

TOTAL  ORGANISMS 

732 

818 

879 

2429 

809.7 

73.9 
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Table  24.  Enumeration  of  macroinvertebrate  taxa  from  Hess  samples  at  the  Missoula  Site,  January  23, 
1997 

Taxon 

Sample  1 

Sample  2 

Sample  3 

SUM 

Mean 

SD 

COLEOPTERA 

0 

0 

3 

3 

1.0 

1.7 

DIPTERA 

40 

87 

110 

237 

79.0 

35.7 

EPHEMEROPTERA 

35 

36 

19 

90 

30.0 

9.5 

PLECOPTERA 

47 

55 

30 

132 

44.0 

12.8 

Claassenia  sabulosa 

10 

6 

5 

21 

7.0 

2.6 

Pteronarcys  californica 

0 

4 

0 

4 

1.3 

2.3 

TRICHOPTERA 

14 

32 

22 

68 

22.7 

9.0 

Hydropsyche  cockerelli 

2 

12 

4 

18 

6.0 

5.3 

Hydropsyche  occidentalis 

8 

16 

10 

34 

11.3 

4.2 

Arctopsyche  grandis 

1 

0 

1 

2 

0.7 

0.6 

TOTAL  ORGANISMS 

136 

210 

184 

530 

176.7 

37.5 

44 


Table  25.  Enumeration  of  macroinvertebrate  taxa  from  Hess  samples  at  the  Missoula  site,  August  12, 
1997 

Taxon 

Sample  1 

Sample  2 

Sample  3 

SUM 

Mean 

SD 

COLEOPTERA 

6 

0 

3 

9 

3.0 

3.0 

DIPTERA 

14 

7 

44 

65 

21.7 

19.7 

EPHEMEROPTERA 

89 

107 

98 

294 

98.0 

9.0 

PLECOPTERA 

22 

16 

16 

54 

18.0 

3.5 

Claassenia  sabulosa 

9 

10 

7 

26 

8.7 

1.5 

Pteronarcys  californica 

11 

6 

7 

24 

8.0 

2.6 

TRICHOPTERA 

268 

276 

233 

777 

259 

22.9 

Hydropsyche  cockerelli 

201 

157 

164 

522 

174 

23.6 

Hydropsyche  occidentalis 

25 

25 

10 

60 

20 

8.7 

Arctopsyche  grandis 

7 

31 

21 

59 

19.7 

12.1 

TOTAL  ORGANISMS 

399 

406 

394 

1199 

399.7 

6.0 
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percentage  was  only  9.7%.  The  1997  year-class  of//,  cockerelli  may  have  been  particularly  strong 
relative  to  other  caddisflies,  but  the  differences  in  density  may  have  also  been  due  to  a  differences 
in  the  emergence  times  of  the  two  previous  year  classes.  The  fact  that  the  H.  cockerelli  were  larger 
in  August  1997  (mean  of  1 1.3  mm  total  length,  minimum  of  10.0  mm)  than  in  July  1996  (mean  10.3 
mm,  minimum  7.5  mm)  suggests  that  the  1996  year-class  had  emerged  and  reproduced  earlier  than 
the  1 995  year-class  and  the  new-generation  larva  were  large  enough  to  be  sampled  or  were  utilizing 
locations  that  we  sampled. 

In  summary,  the  trends  in  seasonal  insect  abundance  seemed  to  follow  the  natural  pattern 
expected  for  temperate-region  streams.  The  one  exception  to  this  was  the  sharp  drop  in  total  insect 
numbers  from  November  1 996  to  January  1 997,  although  we  suspect  that  local  ice  scour  (not  from 
above  Milltown  Dam)  was  responsible. 


Comparison  with  USGS  Data 

There  were  three  instances  in  1996  and  1997  where  the  USGS  sampled  insects  and  sediments 
from  the  same  site  and  at  nearly  the  same  time  as  samples  collected  for  this  study.  In  August  1996, 
the  USGS  sampled  at  the  Alberton  site  about  two  weeks  after  sampling  was  done  for  this  study, 
while  in  1 997,  the  USGS  sampled  both  the  Alberton  and  East  Missoula  sites  about  one  week  earlier 
than  we  did.  The  USGS  analyzed  2-5  replicate  samples  for  each  insect  species  or  sediment,  while 
we  had  2-3  replicate  samples  on  all  occasions  for  sediments  and  Claassenia  and  once  for 
Arctopsyche.  However,  we  only  analyzed  single  samples  for  Hydropsyche  cockerelli  and  in  2  of 
3  instances  for  Arctopsyche,  and  these  single  samples  preclude  reliable  comparisons  between  the 
studies. 

The  one  instance  where  we  analyzed  three  samples  of  Arctopsyche  grandis  was  East  Missoula 
in  1997,  and  the  mean  copper  concentration  in  this  study  was  52.5  ug/g,  within  0.5%  of  the  mean 
value  of  54.6  ug/g  calculated  by  the  USGS  (Table  26).  Mean  values  for  Claassenia  were  always 
somewhat  higher  for  the  USGS,  ranging  from  1%  higher  at  East  Missoula  in  1997  to  30%  higher  at 
Alberton  in  1996.  The  copper  concentrations  calculated  by  the  USGS  for  sediments  were  always 
higher  than  those  in  this  study.  The  mean  USGS  value  was  25%  higher  at  Alberton  in  1996  (136 
vs  108.6  ug/g),  23%  higher  at  Alberton  in  1997  (192  vs  156  ug/g),  and  214%  higher  at  East 
Missoula  in  1997  (516  vs  241  ug/g). 

With  the  exception  of  the  sediment  samples  at  East  Missoula  in  1997,  the  copper  concentrations 
calculated  for  this  study  are  in  good  agreement  with  those  generated  by  the  USGS.  Even  differences 
up  to  30%o  are  probably  not  significant,  given  the  standard  deviation  values  associated  with  each 
mean.  In  addition,  some  of  the  differences  may  be  due  to  variability  in  laboratory  preparation  and 
analysis.  In  this  study,  blind  duplicates  submitted  to  the  laboratory  varied  by  as  much  as  19%. 
Finally,  differences  in  the  exact  locations  of  sampling  may  also  have  contributed  to  the  observed 
variations  in  copper  levels.  The  sediments  we  sampled  at  Alberton  in  1997  were  from  a  spot  several 
miles  upstream  from  where  the  USGS  sampled.    Also,  some  of  the  insects  that  we  sampled  at 
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Alberton  in  1997  were  from  the  same  riffle  sampled  by  the  USGS,  but  some  were  from  riffles  up  to 
a  mile  away. 

Even  with  the  differences  between  the  studies,  the  USGS  data  are  consistent  with  the  finding  in 
this  study  that  copper  concentrations  in  sediments  increased  with  time  at  Alberton.  Between  August 
1996  and  August  1997,  the  USGS  showed  copper  in  sediment  to  increase  41%,  from  136  ug/g  to  192 
ug/g;  likewise,  in  this  study  the  sediment  copper  increased  44%  over  roughly  the  same  time  period. 
The  USGS  samples  at  Turah  showed  that  this  phenomenon  was  not  restricted  to  the  river  below 
Milltown  Dam.  Between  August  1 996  and  1 997,  the  sediment  copper  at  Turah  increased  78%,  from 
356  ug/g  to  635  ug/g. 
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Table  26.  Comparison  of  mean  (+  SD)  copper  concentrations  (ug/g  dry  weight)  in  aquatic  insects 

and  fine-grained  sediments  in  August  1996  and  1997  as  measured  by  the  USGS  (Dodge  et 
al.  1997,  Dodge  et  al  1998,  Hornberger,  pers.  comm.)  and  in  July  1996  and  August  1997 
from  this  study.  Sample  size  is  in  parenthesis. 

Site/                                    Hydropsyche          Arctopsyche          Claassenia 

Collector                            cockerelli                grandis                sabulosa               Sediment 

1996 

Turah  (USGS) 

50.8±11.2(3) 

40.2±6.6  (5) 

45.3±7.6  (3) 

356±18(3) 

E.  Missoula  (This  study) 

39.7(1) 

19.8(1) 

51.2±17.2(3) 

216.7±103(3) 

Missoula  (This  study) 

48.6(1) 

10.7±2.6(2) 

42.2±10.0(3) 

189.7±34(3) 

Aiberton  (TMs  study) 

28,5(1) 

20.1(1) 

43.1+4.2(3) 

108.6+13(3) 

Aiberton  (USGS) 

y 24.5+2.8  (3) 

13.1+2.7(3) 

56.1+8.3(3) 

136+6(3) 

1997 

Turah  (USGS) 

109±7(3) 

92.9±25  (5) 

64.2±18.6(3) 

635±31  (3) 

E.  Missoula  {This  study)   ' 

67.£(1) 

52:5+10.5(3) 

43.8+5.3  (3) 

241.0+42(3) 

E.  Missoula  (USGS) 

89.2+6.4(3) 

54,6+20(3) 

44.4+9.3  (2) 

516+62(3) 

Missoula  (This  study) 

46.3  (1) 

24.7±1.0(2) 

42.0+14.8  (3) 

335.3±51  (3) 

Aiberton  (This  study) 

37.4(1) 

16.4(1) 

44.1+2.8(2) 

156.0+3(3) 

Aiberton  (USGS) 

411+4.7(3) 

28.6+6.9(4) 

49.7+7.7(2) 

192+17(3) 
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